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Abstract 

In this work, we have considered variable G in flat 
FRW universe filled with the mixture of dark energy, 
dark matter and radiation. If there is no interac- 
tion between the three fluids, the deceleration parame- 
ter and statefinder parameters have been calculated in 
terms of dimensionless density parameters which can 
be fixed by observational data. Also the interaction be- 
tween three fluids has been analyzed due to constant 
G. The statefinder parameters also calculated in two 
cases: pressure is constant and pressure is variable. 
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1 Introduction 

The present acceleration of the universe as favored 
by the Supernovae type la data can be explained 
by some exotic matter dominated in the present 
universe which violates the strong energy c o nditio n 
is termed as da rk energy iPerlmutter et all (|l998l ): 
Riess et all (119981 ). This dark energy has the prop- 
erty that it has positive e nergy a nd sufficient negative 



R. R. Caldwelll (I2002D : IR. R. Caldwell et al 



pressu re 

( 2OO2I ). Dark energy occupies about 73% of the en- 
ergy of our Universe, while dark matter about 23% 
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and the usual baryonic matter 4%. There are differ- 
ent candidates obe y the propert y of da r k energy given 



_ quintessence 
k-essence 



Peebles et al (2002D: CaldweU et al 



Armendariz et al| (2000), tachyon 

IA. SenI (l2002l). phan tom Ir,. R. Caldwelll (l200^ . ghos t 
condensate lArkani et a i (|2004h: If iazza et a| (j20Q41 ). 



quintom iFeng et all (120051): iGuo et all (j2005|) . brane 
world models Sahni. Shtaiioy J 20031) and Chaplygin gas 



models iKamenshchik et all (|200ll ). 

Einstein's field equations have two parameters — 
the Newton's gravitational constant and the cosmo- 
logical constant. The Newton's gravitational constant 
G plays the role of a coupling between geometry and 
matter in the Einstein field equations. In an evolv- 
ing universe, it appears natural to loo k at this "con- 
stant" as a function of time. Dirac F. A. M. Dirac 



( 19371 ) proposed for the first time the idea of a vari- 
able G on certain physical grounds. He has shown 
that G oc but his model ran in some difScu l- 

ties. Some aut hors IA-M. M. Abdel RahamanI |l990t ): 
C. Masj |l994l ) have shown that G is an increas- 
ing function of time. Many other extensions of 
Einstein's theory with time dependent G have also 
been proposed in order to achieve a possible unifi- 
cation of gravitation and elementary particle physics 
or to incorporate Mach's princ ip le in general relatiy - 
itv iHovle. Narlikaj (Il964l. Il97ll): iBrans. Dick3 (Il96ll). 



Canuto and Narlikar ICanuto. Narlikarl (jl980l ) have 
shown that the G-varying cosmology is consistent with 
whatever cosmological observations presently available. 
According to Di rac's large numb e rs hyp othesis, G/G ^ 
Hubble rate H Iv. N. Melnikovl |2009l ). Observations 
of Hulse- Taylor binary pulsar B 19 13-1- 16 gives the esti- 
mate < G/G ~ 2 ± 4 X 10-i2yr 



C. S. B. KogI3 



( 20061 ) and helioseismological data gives the boun d 



< G/G ~ 1.6 X IQ-i^ r -i ID B. Guenthed (|199 



S. Rav. U. Mukhopadhvavl |2007t ). Several works 



variable G have been studied in last few decades 
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K.D. Krori et all ( 2000 ): 



1998 



2001 



2002 



20081) . 



M. Jamil U. DebnathI (l201lh : 
M. Jamill |2010f ): lA. I. ArbabI i 

In this work, we have considered the universe is filled 
with radiation, dark matter and dark energy with and 
without interactions. We now follow the method of the 



ref lM. Jamill (|2010[ ) in three fluid system. The dimen- 
sionless density parameters and statcfindcr parameters 
have been calculated due to variable G in three fluid 
system of flat FRW universe. 



pressure i.e., Pm = 0. Also px and p,. are the pres- 
sures of dark energy and radiation respectively. Now 
the equation of state for radiation is given as = \Pr- 
Now assume that the equation of state for dark energy 
is Px = wpx- Next we study the non-interacting and 
interacting situations. 



3 Non-Interacting case 



2 Three Fluids with Varying Gravitational 
Constant 

The isotropic, homogeneous and flat FRW model of the 
universe is described by the line element 

ds^ = dt^ - a^{t) [dr^ + r^{d0^ + sin^ Odc/)^)] , (1) 

where a(t) is the scale factor. The energy-momentum 
tensor is given by 



T^,i, = {p + p)Uf,u„ -pgf^i,, 



(2) 



where is the four velocity satisfying u'^u^ = 1. Here 
p and p are respectively the energy density and pressure 
of the fluid in the universe. 
The Einstein's field equations are given by 



(3) 



where Rf^,, 5^,/ and R are Ricci tensor, metric tensor 
and Ricci scalar respectively. Here we consider gravita- 
tional constant G as a function of cosmic time t. Now 
assume that the universe is filled with the mixture of 
three fluids (radiation, dark matter and dark energy), 
so from the equations (1), (2) and (3) we have the Ein- 
stein's field equations as 



and 



-P, 



(4) 



H = -AnGit){p + p). (5) 
Also the conservation equation is given by 



3H{p + p)=0, 



(6) 



where, p ^ Px + Pm + Pr and p = Px + Pm + Pr are 
the total energy density and pressure of the fluid in 
the universe. Here, PxiPm and p^ are the energy densi- 
ties of the dark energy, dark matter and radiation flu- 
ids respectively. For dark matter, we choose negligible 



If there is no interaction between dark energy, dark 
matter and radiation then from equation (6), we can 
write the continuity equations for dark matter, dark 
energy and radiation as 



Pm + 3Hp,n = 0, 
Px + 3ff (1 + w) Px = 0, 
p^ + 4Hpr = 0. 



(7) 
(8) 
(9) 



Now we define the dimensionless density parameters in 
the form 



Px 



n^^f^ (10) 



, J ^2; , ^^m — ; ^ 'T — 

Per Per Per Per 

where pcr is the critical density, so we obtain 

n = nx + ^im + ^r- (11) 

The deceleration parameter q ^ — can be expressed 
as in terms of density parameters in the following 



From above we obtain 



(12) 



q - [{I + 3w) fix + Swflx + flm + ^r] ■ (13) 



From (9), we have 
_P P_- 

i ' 2 Pcr^ 

Per Per 

where 

f^H g\ 

Per-Per (^^^ " G j ' 

Furthermore, 

H=-H'^{l + q). 
So from (14) we have 
p,,^-Hp,r{2{l + q) + AG), 



(14) 



(15) 



(16) 



(17) 
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where Aq = G'/G, G = HG'. Now from (13) we get 



Pc 



QH [2 (! + (?) + AG] 



From equations (7) to (10), we yield 
= n„M [-1 + 2q + AG] , 

= n^H [-1 - 3w + 2q + AG] , 
and 



QrH [-2 + 2q + AG] . 



(18) 

(19) 
(20) 
(21) 



Using (12) and (22), we can write equation (25) in the 
form 

r= ^[{i + 3w)n, + a,n + nrf 

(-1 + 2(7 + AG) 



(26) 



Substituting equations (19) - (21) in equation (13) and 
after simpUfying, we get 

1 f 3 

q = -n,nH (-1 + 2q + AG) + n^H I (1 + 3w) + -w 

(-1 - 3w + 2q + AG) + QrH (-2 + 2q + AG)] (22) 

We also determine the dimensionless pair of cosmo- 
logical diagnostic pair | r, s} dubbed as sta tefinder pa- 
rameters introduced by ISahni et~ai (2003). The two 
parameters have a great geometrical significance since 
they are derived from the cosmic scale factor alone. 
Also this pair generalizes the well known geometrical 
parameters like the Hubble parameter and the deceler- 
ation parameter. The parameter r forms the next step 
in the hierarchy of geometrical cosmological parameters 
H and q. 

The diagnostic pai r has been u sed fo r holographic 
dark e nergy model Seta,re et all ( 2007 ); Jamil et"al 



|(1 + 3w) + -wj X 
(-1 -3w + 2q + AG) - nr (-2 + 2q + AG)] 
Substituting this in equation (24), we get 

s = { [(1 + 3w) + a,n + 

+ [(-1 + 3w) + O,] - 0,„ {2q + AG) 



-2n^ |(1 3w) + -wj X (-1 - 3w + 2q + AG) 
~2nr i~l + 2q + AG)]} 



(27) 



2009 1) JSetar e fc Vaeenasl (l2009h : 



Setare (200e 

201 Ih , Chaplygin gas lJamil fc DebnathI (|2011[ ) 
201 ll i. 



4 Interacting Case 

Interacting models where the dark energy weakly in- 
teracts with the dark matter have also been studied 
to explain the evolution of the Universe. This models 
describe an energy flow between the components. To 
obtain a suitable evolution of the Universe an interac- 
tinn is often assumed such that the decay rate should 



The parameters are given by 



and 



1 



3(g-i/2) 



For varying gravitation al constant G, t he definition of 
s can be generalized to M. Jamil ( 2010l) 



Setare fc JaaiiA^p^,-f;^noi iq the present value of the Hubble pa- 
Debnath fcr Jjgglll- for good fit to the expansion history of the 
Univ erse as determined b y the Supernovae and CMB 
data Im. S. Bergej |2006l) . These kind of models de- 
scribe an energy flow between the components so that 
no components are conserved separately. First, we as- 
sume that the dark matter component is interacting 
with dark energy component, so the continuity equa- 
tions of dark matter and dark energy are 



(23) 



(24) 



p„i + 'iHp„i = Q, 
3H{l + w) p, = -Q', 



(28) 
(29) 



From (23) , we can write r in terms of H and q as 



2q' 



q 



(25) 



where Q and Q' in order to include the scenario in which 
the mutual interaction between the two principal com- 
ponents of the universe leads to some loss in other forms 
of cosmic constituents. In this case, we have assumed 
Q ^ Q', so from (6), we have the continuity equation 
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for radiation fluid aslCruz et all (12008) : Dcbnathl (|2010|): 
Jamil et ail |2010l ): Ijamfl fc RahamaiJ |2009l) 



AHpr = Q' -Q. 



(30) 



If Q < Q', then part of the dark energy density decays 
into dark matter and the rest in the radiation fluid. 
But if Q > Q', then dark matter receives energy from 
dark energy and from radiation. We are taking about 
in this case that dark energy decay into dark matter 
(or vice versa, depending on the sign of Q) and ra- 
diation. Assume, the interaction terms Q and Q' are 
Q = Sllii? and Q' = 3112 which measure the strength 
of interactions where Hi and 112 have the dimension of 
density. Now assume that G is constant. So Differenti- 
ating equation (5), we have 

H = -At:G {p + p) . (31) 

Also the deceleration parameter q can be expressed as 

1 3 p , , 

1^-7;- t;- (32) 

2 I p 

• Case-I: when p = 0. Equation (31) reduces to 



H = -AttGH 



- 3pm - 3(1 - W - W ]p: 



^' A 
-W — Px - 'ipr 

Id 



Dividing by we yield 

_H_ _ {{w + wIl2+W^)px + {-l)pr} 

m p 



Now equation (25) can be written as 
2. 



H 



Using (34), the above equation becomes 

25 9p (w + wU2 + up) .px 

= ^ + T" - ^- 

2 4p p 

Also from (23), we obtain the expression for s as 

23 9p-36(i«-(-uin2+to^)pa; 
_ ~ ' 4p 



(33) 



(34) 



(35) 



(36) 



(37) 



• Case-II: when p ^ then for p — Px ^ Pr — 

wpx + \pr we get (from (31)) 

^ = -47rG |-3pm + ^-3 ~&w ~ 3w^ + j^"^ Px 
16 



(38) 



p^-(3w-i)n2-ni 



(1 - w) - W — px - 4pr 



Dividing by we obtain 
H 3 1 



m 2p 
16 



-3/5m + [ -3 — Gw — 3w'^ + —w ) p, 



— - (3u. - 1) n2 - Hi X 

Q' 

(1 - W) - W — px - 4pr 



(39) 



In this case, the expressions of r and s become 

7 Qpm ^(3 9 \px 13 pr 3D2 , , 

r = - - - — + -Di + 7:^] — — , 40) 

2 2 p \2 2 J p 2 p 2 p 



and 



3f-l-fH 



(41) 



where 



Di = —3 — &w — 3w -I- —w 
H 

D2-- (3^-1)02 -Hi. 



5 Discussions 

In this work, we have considered variable G in flat FRW 
universe filled with the mixture of dark energy, dark 
matter and radiation. If there is no interaction be- 
tween the three fluids, the deceleration parameter and 
statefinder parameters have been calculated in terms 
of dimensionless density parameters which can be fixed 
by observational data. Also the interaction between 
three fluids has been analyzed due to constant G. If 
Q < Q\ then part of the dark energy density decays 
into dark matter and the rest in the radiation fluid. 
But if Q > Q' , then dark matter receives energy from 
dark energy and from radiation. We are taking about 
in this case that dark energy decay into dark matter (or 
vice versa, depending on the sign of Q) and radiation. 
The stateflnder parameters also calculated in two cases: 
pressure is constant and pressure is variable. In the lit- 
erature, the diagnostic pair has been calculated for the 
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dark energy interacting with dark matter iJamil et a] 
( 2010t ). The interaction of these two dark fluids with 
a third component has been studied in hterature, and 
hence motivated by this, we calculated the statcfinder 
parameters for the triple fluid interacting model. This 
model will be useful for investiga t ing fu rther the triple 
coincidence problem I Arkani et all (|2004l ). 
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